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Abstract 

We study the 0[a t a s + af ) two-loop corrections to the minimization conditions of the MSSM 
effective potential, providing compact analytical formulae for the Higgs tadpoles. We connect 
these results with the renormalization group running of the MSSM parameters from the grand 
unification scale down to the weak scale, and discuss the corrections to the Higgs mixing 
parameter /j, and to the running CP-odd Higgs mass m A in various scenarios of gravity- 
mediated SUSY breaking. We find that the 0(ceta s ) and O(a^) contributions partially cancel 
each other in the minimization conditions. In comparison with the full one-loop corrections, 
the 0(ata s + af) two-loop corrections significantly weaken the dependence of the parameters 
/i and m A on the renormalization scale at which the effective potential is minimized. The 
residual two-loop and higher-order corrections to fj, and m A are estimated to be at most 1% 
in the considered scenarios. 
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1 Introduction 



One of the most attractive features of the Minimal Supersymmetric extension of the Standard 
Model (MSSM) PQ, is the fact that it provides a mechanism for breaking radiatively the elec- 
troweak gauge SU(2)l x U(1)y symmetry down to £/(1)em- It was first shown [2] that a 
supersymmetry (SUSY) breaking term for the gluino can induce an effective potential which 
spontaneously breaks the electroweak symmetry. At the same time, a mechanism relying on 
the renormalization group evolution from a grand unification (GUT) scale M GUT down to the 
weak scale was proposed .3 . In this framework, at the scale M GUT , the parameters entering 
the scalar potential of the MSSM obey simple boundary conditions dictated by the underlying 
theory of SUSY breaking, and the electroweak symmetry is unbroken. When the parameters 
are evolved down to the weak scale by means of the MSSM renormalization group equations 
(RGE), which amounts to resumming the leading logarithmic corrections to all orders, the soft 
SUSY-breaking mass rn 2 H2 is driven towards negative values, due to corrections controlled by 
the top Yukawa coupling ht- This helps to destabilize the origin in field space, so that the Higgs 
fields acquire non-vanishing vacuum expectation values (VEVs) and the electroweak symmetry 
is spontaneously broken. Although the studies in Refs. [3 IS] where differing on the initial bound- 
ary conditions, the result was one: the radiative electroweak symmetry breaking (REWSB) takes 
place if the top Yukawa coupling is large, such that 60 GeV < m t < 200 GeV, with the upper 
bound coming from the requirement that ht remains in the perturbative range up to the GUT 
scale. It could be a coincidence that the top quark is found at the Tevatron to have mass around 
175 GeV, but certainly this is consistent with the REWSB mechanism in the MSSM. 

In the RGE-improved potential of the MSSM employed at tree level, the VEVs of the Higgs 
fields, and the occurrence of spontaneous symmetry breaking itself, depend critically on the renor- 
malization scale at which the parameters entering the potential are computed; an inappropriate 
choice of that scale can lead to results that are even qualitatively wrong. In fact, the electroweak 
symmetry is either broken or unbroken, independently of the renormalization scale choice, and 
the critical behavior described above is just an artifact of the tree-level approximation. The 
correct way of determining the ground state of the theory is to minimize the Coleman- Weinberg 
effective potential [I], i.e. the tree-level potential plus a correction coming from the sum of all 
the one-loop diagrams with zero-momentum external lines. Since Refs. this procedure has 
become standard in the renormalization group analyses of the MSSM (for early examples see 

Refs. mini). 

The effective potential is also a useful tool for computing the leading corrections to the 
MSSM Higgs masses, both at the one loop and the two loop QUI QU 113 H3 HI level 1 , in 
the approximation of zero external momentum. The leading one-loop corrections are 0{at), i.e. 
they are controlled by the top Yukawa coupling h t = \f^naf For stop masses of 0(1 TeV), such 
corrections increase by 40-60 GeV the mass of the lightest Higgs boson (which at tree level 
must be lighter than m z ), allowing it to escape the direct searches at LEP. Also, the leading two- 
loop corrections have sizeable effects: the 0{ata s ) corrections, controlled by the strong gauge 
coupling g s = ^/A^ra^ , typically reduce by 15-20 GeV, whereas the 0(pq) ones may increase 
it by up to 7-8 GeV. 

1 Other two-loop computations of the MSSM Higgs masses have been performed in the renormalization 
group 15 and diagrammatic |16| approaches. 
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Motivated by the relevance of the 0{oct(x s + af ) two-loop corrections in the case of the 
Higgs masses, we study in this paper the effect of the same corrections on the electroweak 
symmetry breaking conditions. The contributions to the two-loop MSSM effective potential 
that are relevant to the (D(ata s + af) corrections have been discussed in Refs. IT2*] . 
and a complete computation of the two-loop effective potential has been presented in Ref. |17| . 
However, practical studies of REWSB usually require explicit formulae for the Higgs tadpole 
diagrams, i.e. the first derivatives of the effective potential with respect to the Higgs fields. Such 
formulae are presently available at the one-loop order [7j, but they have not been presented 
so far at the two-loop order. Using the techniques developed in Ref. [T2|, we compute in this 
paper explicit and compact analytical expressions for the two-loop 0(ata s + af) part of the 
tadpoles. As a byproduct from our 0(ata s ) corrections we obtain also the 0(aba s ) corrections, 
that are relevant for large values of tan/3. Once we assume that the electroweak symmetry is 
indeed broken, giving rise to the observed value of the Z boson mass, the corrections to the 
tadpoles translate into 0(ata s + af) corrections to the values of fi, the Higgs mass term in the 
superpotential, and m A , the running mass of the A boson. We discuss the effect of our two-loop 
corrections in the framework of gravity mediated SUSY breaking [1B|. also denoted as minimal 
supergravity (mSUGRA), referring in particular to various "benchmark" scenarios suggested at 
Snowmass |19j . We find that the inclusion of the 0{ata s + af) corrections significantly improves 
the renormalization scale dependence of the results, and that partial cancellations occur between 
the 0(ata s ) corrections and the 0(af) ones. Our corrections are also required for consistency in 
the 0(a t a s + af) two-loop computation of the MSSM Higgs masses, if the input parameters are 
computed via renormalization group evolution from a set of high energy boundary conditions. 

The paper is organized as follows: in section 2 we recall the basic concepts of radiative 
electroweak symmetry breaking, and introduce some notation which will be used in the rest 
of the paper; in section 3 we describe the main features of our 0(ata s + af) computation of 
the two-loop tadpoles; in section 4 we discuss the numerical effect of our corrections, and we 
show how they improve the dependence of [i and m A on the renormalization scale at which the 
effective potential is minimized; section 5 contains our conclusions. In addition, we present in the 
appendix A some useful formulae for the integrals entering the two-loop effective potential, and 
in the appendix B the explicit analytical formulae for the 0(ata s ) part of the corrections. The 
formulae for the O(af) part are indeed rather long, thus we make them available, upon request, 
in the form of a computer code 2 . 



2 Radiative electroweak symmetry breaking 

We start our discussion from the tree-level scalar potential of the MSSM, that reads, keeping 
only the dependence on the neutral Higgs fields H® and H®: 



V = A + m 2 



Hi 



Hi 



+ mi [H^Hl + h.c.) + 



9 2 + 9' 2 



(\Hl\ 



(1) 



where: A is a field-independent vacuum energy; m 2 = m 2 Hi +ii 2 , m 2 = m 2 H2 +n 2 (we assume /i to 
be real, neglecting all possible CP-violating phases); m 2 H , m^ 2 and m 2 are soft SUSY-breaking 

2 E-mail: slavich@mppmu.mpg.de 
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masses; g and g' are the SU(2)l and U(1)y gauge couplings, respectively. At the classical level, 
the mass parameters entering Vq must satisfy the following conditions: 



mi+m 2 > z|m 3 | , rrij m 2 < m 3 . (zj 

The first condition guarantees that the potential is bounded from below; the second condition 
destabilizes the origin in field space, making sure that the neutral components of the Higgs 
fields acquire non- vanishing VEVs (H®) = vi/y/2 and (H§) = v 2 /V2. It is not restrictive to 
choose m 3 real and negative 3 , so that v± and v 2 are real and positive, and the neutral Higgs 
fields can be decomposed into their VEVs plus their CP-even and CP-odd fluctuations as Hf = 
( Vi + Si + iPi)/V2. 

Since the parameters entering Vq are taken as "running" ones (i.e., they vary with the renor- 
malization scale), also the validity of the conditions in Eq. © depends on the scale, as well as 
the numerical values of v\ and v 2 . As discussed in Ref. jS], the minimization of the tree-level 
potential may lead to grossly inaccurate results, unless the renormalization scale is chosen in 
such a way that the radiative corrections to the scalar potential are small. To obtain the correct 
results, one should rather minimize the effective potential V e s, defined as: 

V cS = V + AV, (3) 

where AV contains the radiative corrections to the scalar potential Vo- The minimization con- 
ditions for V e g can be written as: 



1 dVeB 



v x dSi 
1 dV cS 



t>2 dS 2 

where the "tadpoles" Si and T, 2 are defined as: 

_ 1 dAV 

Zj,; = — 



2 / 2 

2 , 2,5+5' / 2 2\, 2 W 2.^ n /.x 

m H 1 +A t + 1 ( v i -v 2 ) +m 3 — + Si =0, (4) 

2 / 2 

2 i 2 i 9 + 9 I 2 2\ i 2 v l i n fz\ 

m H 2 +V + iv 2 -v 1 )+m 3 — + S 2 =0, 5) 

4 v 2 



Vi dSi 



(6) 



In principle, a renormalization group study of the MSSM should start from some large scale 
M GUT , where the input parameters have a simple structure dictated by the underlying theory of 
SUSY breaking, and the electroweak symmetry is unbroken. The parameters are then evolved, 
by means of appropriate renormalization group equations, down to some lower scale, where the 
electroweak symmetry breaking occurs and the VEVs v\ and v 2 can be obtained by solving 
Eqs. A set of high-energy input parameters is then acceptable if it leads to the correct 

value of the squared running mass for the Z boson, m| = (g 2 + g' 2 ) (v \ + v 2 )/4. However, in most 
practical applications of the renormalization group procedure, it is more convenient to assume 
that there is successful electroweak symmetry breaking, and trade two of the high-energy input 
parameters for v\ and v 2 (or, equivalently, for v 2 = v 2 + v 2 and tan/3 = V2/V1). Eqs. ©-© can 

3 Our conventions differ by a sign in the parameters fj, and m| with respect to those used in the second paper 
of Ref. [T]. 
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thus be rephrased into the following conditions among the parameters at the weak scale: 



m| m 2 Hl + Si - (m 2 H2 + S 2 ) tan 2 (3 
2 + tan 2 /?-l ' 1 j 

X - sin 2/? (m 2 Hl + m 2 H2 + 2 fi 2 + Si + S 2 ) , (8) 

i.e., the terms proportional to Si and S 2 in the above equations can be viewed as the radiative 
corrections to the values of /x 2 and m\ obtained from the requirement of successful electroweak 
symmetry breaking. We recall that m 2 is related to the squared running mass for the A boson 
through m 2 = —2 m 2 ;/ sin 2(3. Notice also that the sign of [i is not fixed by Eq. (J7J), and it 
must be supplemented as an additional input quantity. If the right side of Eq. ("""} is such 
that jj? is negative, then our choice of input parameters is inconsistent, and the electroweak 
symmetry fails to be broken. We remark in passing that the choice of the input parameters 
is constrained by further requirements: it must lead to a spectrum of physical masses for the 
MSSM superpartners and Higgs bosons compatible with the present experimental lower bounds, 
and such that the lightest supersymmetric particle (LSP) is electrically neutral; it must satisfy 
phenomenological constraints coming from radiative B-meson decays, muon anomalous magnetic 
moment and cosmological relic density; finally, it must guarantee that the MSSM scalar potential 
is bounded from below and does not lead to charge and color breaking minima. However, a 
detailed study of the (theoretically and experimentally) allowed regions in the MSSM parameter 
space goes beyond the scope of this paper, and will not be pursued in the following. 

The full one-loop corrections to the REWSB conditions have been extensively discussed in the 
literature [7] in the framework of the mSUGRA scenario. The dominant one-loop contributions 
to Si and S 2 come to a large extent from the top/stop (and, for large tan/3, bottom/sbottom 
and tau/stau) diagrams. The contributions of the diagrams involving charginos and neutralinos 
can also be sizeable and comparable to the top/stop ones in some regions of the parameter space, 
while the Higgs and gauge bosons and the first two generations of (s)quarks and (s)leptons give 
only subdominant corrections. 

In the following sections we provide explicit analytical formulae for the two-loop top/stop 
contributions to the tadpoles Si and S 2 , resulting into 0(otta s + a 2 ) corrections to /u 2 and to the 
running mass m 2 . In analogy with the case of the Higgs masses, we expect such corrections to be 
the leading ones at the two-loop level, giving rise to sizeable effects at least in some regions of the 
MSSM parameter space. In addition, most public codes that compute the MSSM mass spectrum 
from a set of unified parameters at the scale M GUT , such as SuSpect [20], SoftSusy ["""Q) SPheno 
[2*2*] and FeynSSG [23]) include a two-loop (D(ata s + a 2 ) computation of the Higgs masses, 4 but 
employ one-loop results for the tadpoles (see also Ref. ["""5] for a recent discussion). Since m 2 
enters the tree-level mass matrix of the CP-even Higgs bosons, the 0(ata s + a 2 ) corrections to 
m\ should be included in those codes for consistency. 

4 Another widely used public code, Isajet 7.58 12 4| . relies on a one-loop effective potential computation of the 
Higgs masses. 
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3 Computation of the 0(ata s + af) corrections 

We shall now describe our two-loop, 0(ata s + a 2 ) computation of the tadpoles Si and £2, 
involving the first derivatives of AV with respect to the CP-even parts of the neutral Higgs 
fields [see Eq. 1(5])]. The computation is consistently performed by setting to zero all the gauge 
couplings but g s and by keeping lit as the only non-vanishing Yukawa coupling (with a slight 
abuse of language, in the following we will refer to this approximation as to the gaugeless limit). 
In this limit, the tree-level (field-dependent) spectrum of the MSSM simplifies considerably: 
gauginos and Higgsinos do not mix; the charged and neutral Higgsinos combine into Dirac spinors 
hP and hp 1 with degenerate mass the gaugino masses coincide with the soft SUSY-breaking 
parameters M4 (A = 1,2,3) (among them, only the gluino mass m g = M3 is relevant to our 
calculation); the only massive Standard Model (SM) fermion is the top quark; all other fermions 
and gauge bosons have vanishing masses; besides the top squarks, the only sfermion with non- 
vanishing couplings is the bottom squark bi', the lighter CP-even Higgs boson, h, is massless, and 
the same is true for the Goldstone bosons G and G ± ; all the remaining Higgs states, (H, A, H^), 
have degenerate mass eigenvalues m 2 A . The tree-level mixing angle in the CP-even sector is just 
a = 0- tt/2. 

To begin with, we address the renormalization of the effective potential. In the loop expan- 
sion, the correction to the effective potential can be decomposed as AV = Vu + V21 H j where 

the ellipsis stand for higher loops. Using the Landau gauge and dimensional reduction [^Hl HZj 
in d = 4 — 2e dimensions, and including for later convenience also terms that vanish when e — > 0, 
the unrenormalized one-loop effective potential reads: 



1 3 , M 2 7 3, M 2 1, 2 M 2 tt 2 

— I In — - + e In — - + - In 2 — H 

e 2 Q 2 U 2 Q 2 2 Q 2 12 



(9) 



where M 2 is the matrix of the field-dependent squared masses, and the supertrace of a generic 
function f(A4 2 ) is defined as a sum over the eigenvalues m 2 : 

Sir f(M 2 ) = £(-l) 2 *(2 Si + 1) f(m 2 ) , (10) 

i 

where Sj is the spin of the corresponding particles. In Eq. 0, Q 2 = 47r / u 2 e -7B , i.e. the finite terms 
that are removed together with 1/e in the modified subtraction schemes have been reabsorbed in 
the renormalization scale (the same convention will be adopted in the following). In the gaugeless 
limit described above, only the top and stop contributions to Vu are relevant, giving rise to 0(at) 
contributions to £1 and £2- 

The Feynman diagrams that contribute to the two-loop effective potential V21 and give rise 
to 0{atOi s ) contributions to £1 and £2 are shown in Fig. ^ while the diagrams relevant to 
the 0(a 2 ) contributions are shown in Fig. [2j The corresponding analytical formulae for V21 
can be found e.g. in the last paper of Ref. 5 . These formulae involve two basic integrals, 
I(m 2 ,1712,771^) and J(m 2 ,771?,), that have been evaluated with different methods in Refs. |28l I29j . 
Explicit expressions for I and J in the formalism of Ref. [201 are presented in the appendix A. 

5 The formulae of Ref. are obtained for vanishing CP-odd fields. While requiring some modifications |12| 
for the computation of the CP-even Higgs boson masses in terms of the physical tua, those formulae can be used 
as they stand for the purposes of the present analysis. 
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To carry out the renormalization of Veg, at the two-loop order and in the DR scheme |26U27| . 
we start from the unrenormalized effective potential, written in terms of generic bare parameters 
Xi. Then, we expand the parameters as Xi = xf R + Sxi, where 5xi are purely divergent quantities, 
so that all the poles in 1/e and 1/e 2 are cancelled. After taking the limit e — » 0, the non vanishing 
part of the renormalized effective potential is: 

V* = V (xf) + V$(xr) + V 2 f + 6 u xi , (11) 

where vjp and V$ denote the finite parts of the one-loop and two-loop effective potential, 
respectively, V^f denotes the terms proportional to e in Eq. ©, and 5 u Xi is the coefficient of 
1/e in the one-loop part of the generic counterterm (notice that we need to compute explicitly 
only the one-loop counterterms for the top and stop masses). In Eq. Vq and V^p are 

expressed in terms of DR-renormalized parameters, while the renormalization of the parameters 
entering the two-loop part is irrelevant, amounting to a higher-order (i.e., three-loop) effect. In 
summary, it is possible to define the renormalized two-loop effective potential as: 

V 2e = V^ + ^6 u x i . (12) 

OXi 

We have checked that V 2 e corresponds to the finite part of the potential obtained by replacing 
the integrals / and J in V 2 £ with the "subtracted" integrals I and J, first introduced in Ref. |28| . 
More precisely, this is true only up to terms that give a null contribution to £1 and £2, unless we 
include in V e g also diagrams that do not depend on the Higgs fields (such as, e.g., the one-loop 
diagram involving gluinos and the two-loop diagram involving gluinos and gluons). 

Compact analytical formulae for the derivatives of the renormalized effective potential in the 
gaugeless limit can be obtained with a procedure similar to that of Ref. ^2j. The relevant field- 



dependent quantities are the top mass m t , the stop masses m| and m~ , and the stop mixing 
angle % (the top and stop phases ip and ip, introduced in J2j to take into account the dependence 
on the CP-odd part of the Higgs fields, do not enter the computation of £1 and £2). At the 
minimum of the effective potential, the parameters in the stop sector are related by: 

• oa 2 ™t {A t + 11 cot/3) 

sm20r= 5 K , (13) 

mf -mf 
ti ti 

where At is the soft SUSY-breaking trilinear coupling of the stops [notice that Eq. (|L3J) defines 
our convention for the sign of fj,] . We will use in the following the shortcuts C20 = cos 26^ and 
s 26 — sm 2% • After a straightforward application of the chain rule for the derivatives of the 
effective potential, we get: 

r^£i = m t \x cot [3 s 2 e F , (14) 
v\^ 2 = m t A t s 2 eF + 2m 2 t G. (15) 

The functions F and G are combinations of the derivatives of AV with respect to the field- 
dependent parameters, computed at the minimum of the effective potential: 

dAV dAV 2c 2e dAV 



dm? <9m? m| — m~ dcoa 
ti t 2 h t 2 za 



(16) 
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G dAV + dAV + dAV (17) 
dm? dm? dm? 

The one— loop parts of F and G, giving rise to 0(at) contributions to Si and S2, are easily 
computed from the derivatives of v[i > . In units of N c /(16tt 2 ), where N c = 3 is a color factor, 
they read: 

(m? \ ( m? \ 

In-^-lj-m? ^In-^-lj , (18) 

G U = <(-^-l)+-|(l^-l)-2m|(ln^|-l). (19) 

Although a naive, brute-force computation of the derivatives of the renormalized two-loop poten- 
tial Vn presents no major conceptual difficulties, the number of terms involved blows up quickly, 
giving rise to very long and complicated analytical expressions. However, in the spirit of Ref. |29j . 
it is possible to obtain recursive relations for the derivatives of the integral I(m 2 , m 2 , m 2 ) with 
respect to the internal masses (see the appendix A), that simplify considerably the results. In 
this way, we obtained compact analytical formulae for the two-loop parts of F and G, giving 
rise to 0(ata s + a 2 ) contributions to Si and S2. As a non trivial check of the correctness of our 
computation, we have verified that the quantities /i 2 and m 2 defined in Eqs. 0-© obey the 
appropriate two-loop RGE [3U], specialized to the gaugeless limit: 

d/j 2 3 a t 2 , a ta s 2 9a t 2 rnn\ 



dlnQ 2 4vr n ir 2 n 16vr 2 



dm? 3 at ( m? \ ata s ( m? . \ 9 a? ( m? ^ 



(21) 



We have also checked explicitly that, in the special supersymmetric limit in which all the soft- 
SUSY breaking parameters as well as /j, are set to zero, such that mf = m^ 2 = mt, the two-loop 
parts of F and G are indeed vanishing. 

For illustrative purposes, we present in the appendix B the explicit formulae for the 0(ata s ) 
corrections, valid for arbitrary values of the relevant MSSM parameters. The corresponding 
formulae for the 0(a 2 ) corrections are indeed rather long, thus we make them available, upon 
request, in the form of a Fortran code. 

The computation described above allows us to obtain also the two-loop 0(aba s ) corrections 6 
induced by the bottom/sbottom sector (hf, = being the bottom Yukawa coupling), that 

can be relevant for large values of tan/3. To this purpose, the substitutions t — ► b, v\ «-> V2 (i.e., 
tan/3 <-> cot (5) and Si <-> S 2 must be performed in the corresponding formulae for the 0(ata s ) 
part of the corrections. The complications relative to the on-shell definition of the sbottom 
parameters, discussed in Ref. ^31 > do not arise in this case since we are working in the DR 
renormalization scheme. However, the tan /3-enhanced threshold corrections |5J to the relation 
between hb and the bottom mass mt, must be resummed to all orders 32 in a redefinition of the 
bottom Yukawa coupling (see e.g. Ref. f° r the details). 

6 In contrast, the 0(a^ + atotb) corrections would require a dedicated computation. 
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4 Numerical results 



In this section we discuss the numerical effect of our 0(ata s + a 2 ) two-loop corrections on the 
minimization conditions of the MSSM effective potential. 

For definiteness, we work in the mSUGRA scenario, in which the MSSM Lagrangian at 
the large scale M GUT contains only five independent mass parameters: a common soft SUSY- 
breaking scalar mass mo, a common soft gaugino mass mi/2, a common soft trilinear term Aq, 
the superpotential Higgs mixing parameter hq and its soft SUSY-breaking counterpart Bq (the 
subscript "0" denotes the fact that the parameters are computed at the boundary scale). The 
soft Higgs mixing parameter B has the dimensions of a squared mass, and is defined in such a 
way that in the low-energy Higgs potential of Eq. Q it coincides with m 2 (to avoid confusion, 
we will refer to B as to m 2 from now on) . As anticipated in section rather than providing 
input values for all the five mass parameters at the GUT scale, we assume that the electroweak 
symmetry is successfully broken at the weak scale, and we trade fiQ and m 2 (M GUT ) for the weak 
scale input parameters v and tan j3. 

Before discussing our results, it is useful to describe in some detail the numerical procedure 
for the renormalization group evolution of the MSSM parameters. We start by defining the 
SU(3) x SU(2)l x U{1)y gauge couplings at the weak scale (which we identify with the pole 
Z boson mass, M z = 91.187 GeV), from the running weak mixing angle s 2 = 0.2315, the 
electromagnetic coupling a EM = 1/127.9 and the strong coupling a s = 0.119. The electroweak 
symmetry breaking parameter, v 2 = v\ + v\ , is defined in terms the muon decay constant 
according to the relation v = (y/2 G^)^ 1 / 2 = 246.218 GeV, and then translated to the DR 
scheme by means of the formulae of Ref. 12 . In addition, tan/3 = 1)2/ v\ is taken as an input 
parameter at the weak scale, allowing us to determine v\ and V2- The Yukawa couplings of the 
light SM fermions are obtained from the corresponding masses at the scale Q = 2 GeV, and then 
evolved up to M z by means of the two-loop SM RGEs jSHj- In the case of the bottom coupling, 
the tan /3-enhanced threshold corrections [^ E2] to the relation between hb and are included 
at the scale M z according to the formulae of Ref. 13 . Finally, the top Yukawa coupling is also 
defined at the scale M z through the relation ht = \^2mt/v2, where mt is the DR mass for the 
top quark, obtained from the pole mass Mt = 174.3 GeV by means of the formulae of Ref. 12 . 

The evolution of the MSSM parameters from the weak scale to the GUT scale, and back, 
is performed by means of the full one-loop MSSM RGEs. However, for consistency with our 
two-loop 0(otta s + a 2 ) analysis of the REWSB conditions, we supplement the RGEs with the 
two-loop strong and top- Yukawa contributions as given in Ref. |30| . We also make use of 
the two-loop RGEs in the Landau gauge for the VEVs, following Refs. [171 I34j . As a first 
step, we evolve the gauge and Yukawa couplings from the scale M z up to the scale where the 
SU(2)l x U(l)y gauge couplings g\ and 52 meet, that we identify with M GUT (we do not force 
the strong gauge coupling 53 to meet g\ and 52)- At this scale, which turns out to be of the 
order of 10 16 GeV, we set the input boundary conditions for the soft SUSY-breaking masses 
mo, m 1 / 2 and Aq. The values of fiQ and m^{M GUT ), to be later determined from the REWSB 
conditions, are provisionally set to zero. At this point, we start an iterative procedure: first we 
run the MSSM parameters from M GUT down to some scale Q m im of the order of the weak scale, 
where the values of /i 2 (Qmin) and m 2 (Q m i n ) are computed through Eqs. (0) and (jHJ), with the 
sign of [i supplied as an extra input parameter. Then we run all the parameters, including [i 
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and 771,3, down to the scale M z , which we regard as the end point of the RGE evolution. At 
this scale, we compute the threshold corrections to the top and bottom Yukawa couplings, using 
the newly obtained values of the relevant MSSM parameters. Finally, we run all the parameters 
back to the scale M GUT , where the resulting values for \x and m| are taken as new guesses for 
the corresponding boundary conditions. We iterate the procedure until convergence is reached, 
i.e. the values of /i 2 ((5min) and m^(Q m - m ) obtained from the RGE evolution of /io and m|(M GC/T ) 
coincide with those obtained from the minimization conditions (JJJ) and (|SJ) . If however /j 2 turns 
out to be negative, then our choice of input parameters is inconsistent (i.e., it does not lead to 
successful REWSB) and must be discarded. 

In order to discuss the effect of our 0{otta s + a 2 ) corrections to the REWSB conditions, we 
show in Figs. OHZlthe values of |/j(M z )| and m A (M z ), the latter obtained through m 2 (M z ) = 
— 2 m 2 (M z )/ sin 2(3 , as functions of the minimization scale Q m \n- Stability of the results with 
respect to moderate changes in Qmin (which should anyway lie in the weak range, i.e. between 
M z and a few TeV) indicates that the higher-order corrections not included in the computation 
of fj, and m A are small and can be safely neglected. We will see that in general the inclusion of 
our 0(otta s + a 2 ) corrections significantly improves the scale-dependence of the results. 

In the choice of the input parameters, we refer to the so-called Snowmass Points which 
represent typical "benchmark" scenarios that are commonly investigated in the phenomenological 
analyses of the mSUGRA parameter space. In particular, Figs. EHU correspond to: 

SPS la : m = 100 GeV, m 1/2 = 250 GeV, A = -100 GeV, tan/3 = 10 , n < , 

SPS 3 : m = 90 GeV, m 1/2 = 400 GeV, A = , tan/3 = 10 , fi < , 

SPS 4 : m = 400 GeV, m 1/2 = 300 GeV, A = , tan f3 = 50 , fj, < , 

SPS 5 : m = 150 GeV, ro 1/2 = 300 GeV, A = -1 TeV , tan^ = 5 , /j < , 



respectively. Notice that our convention for the sign of [i differs from the one in Ref. |19| . where a 
discussion on the characteristics of the various scenarios can be found. A further scenario, denoted 
as "Focus point" (SPS2) and characterized by a common scalar mass, mo = 1450 GeV, much 
larger than the common fermion mass, m.1/2 = 300 GeV, has also been suggested in Ref. |19| . 
However, we found that in this scenario the results for /j and m A , including the qualitative effect 
of the various corrections and the occurrence of REWSB itself, depend dramatically on very 
small adjustments of the input value for the top pole mass (e.g. |/j(M z )| varies roughly between 
400 and 100 GeV if Mt is varied between 174.3 and 175 GeV). The extreme sensitivity of the 
SPS2 scenario on the input top mass has already been discussed in Ref. |25| . Since this scenario 
appears to lead to unstable results, we will not consider it further in this work. However, in 
order to investigate the situation in which the common scalar mass is considerably larger than 
the common fermion mass, we show in Fig. [7|a scenario with mo = 1 TeV and m 1( / 2 = 300 GeV 
(the other parameters being chosen as Aq = , tan/3 = 10 and \x < 0). We have checked that 
this "Large mo" scenario is not unreasonably sensitive to small variations in the input top mass. 

In all the plots of Figs. [3HZ1 the minimization scale varies in the range M z < Q m \ a < 2Q*, 
where Q* = (m§ + 4m 2 y 2 ) 1 / 2 is a scale roughly comparable with the squark masses. The dotted 
curves in the upper and lower panels of each figure represent |/i(M z )| and m A (M z ), respectively, 
as obtained by including in the minimization conditions only the one-loop 0(ctt) top/stop con- 
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tributions to Si and S2; the dashed lines include instead the full one-loop computation of Si 
and £2; the dot-dashed lines include in addition the two-loop 0(atac s ) contributions; finally, 
the solid lines include our full two-loop result, i.e. the 0(ata s + a^) contributions to Si and £2. 
In Fig. [2J corresponding to the "large tan/3" (SPS4) scenario with tan/3 = 50, we show also the 
effect of the 0(aba s ) corrections, obtained from the 0{ata s ) ones as described at the end of the 
previous section. The dot-dot-dashed lines in Fig. [5] include the 0(ata s + ai,a s ) contributions, 
while the solid lines represent the full 0(ata s + a^cts + a$) result. The effect of the 0(ctba s ) 
corrections is indeed negligible in the other SPS scenarios, where tan (3 takes on more moderate 
values. In any case we include the 0{a\,a s ) corrections in all the scenarios we investigate here. 

We see from Figs. OHZIthat the inclusion of the two-loop 0(ata s + a^) corrections improves 
the dependence of |/i(M z )| and m A (M z ) on the minimization scale with respect to the full one- 
loop result. The effect is particularly manifest in the case of fi, where, in all the scenarios, 
the two-loop corrected result appears to depend only very weakly (within 3 GeV at most) on 
Qrnin, while the one-loop result tends to decrease for both small and large Qmin- In the case of 
m A , the improvement is less striking: although the two-loop corrected result has in general a 
better scale dependence than its one-loop counterpart, especially for increasing values of Q m i n , 
a small residual scale dependence is visible in most plots when Q m \ n gets close to M z . In any 
scenario the residual uncertainty on m A is never larger than 10 GeV (the latter case occurring 
in the "large tan/3" scenario, where m A {M z ) is around 1200 GeV) and might be due to the 
corrections that we neglect in our two-loop computation of the tadpoles, i.e. those controlled 
by the electroweak gauge couplings and, for large tan/3, those of 0(a^ + a t a.b) ■ From the 
small residual scale dependence visible in Figs. EMU we estimate that the effect of the neglected 
two-loop and higher-order corrections on the parameters /i and m A should be at most of 1%. 

Other interesting observations can be drawn from Figs. I2H3 first of all, at the one-loop 
level, the inclusion of the top/stop contributions only is in general not a good approximation 
of the full result (this has been already observed e.g. in Ref. |35|V Moreover, at the two-loop 
level, a significant compensation occurs between the 0(ata s ) and the 0(a^) contributions to the 
tadpoles, thus including only the former may lead to rather inaccurate predictions (as shown by 
the dot-dashed curves). This partial cancellation between the 0{ata s ) and O(a^) corrections 
is similar to the one occurring in the case of the Higgs masses |10| lll| I12j . Finally, it is worth 
noticing that in the "large tan /?" scenario of Fig.Elthe inclusion of the 0(abCt s ) corrections to the 
tadpoles has a sizeable effect on the minimization scale dependence of m A (M z ) (see the difference 
between the dot-dashed and dot-dot-dashed curves), while it does not affect significantly that 
of HMz)\- 

It is clear from the above discussion that the numerical effect on |^(M Z )| and m A {M z ) of 
the two-loop corrections to Si and S2 depends critically on the choice of the minimization 
scale. In all the plots we find a range of values of Q m in> usually in the vicinity of Q*, for 
which the one-loop and two-loop curves are close to each other, implying that the effect of the 
0(ata s + al ) corrections is small. On the other hand, for values of Q m i n far from this optimal 
choice, the omission of the two-loop corrections can lead to an error of several (possibly, tenths 
of) GeV, especially in the case of /i. Thus, the proper inclusion the 0(ctta s + ) two-loop 
tadpole corrections on the top of the full one-loop ones allows us to obtain more precise and 
reliable results for /x and m A , i.e. results that do not depend on a preconceived choice of the 
minimization scale for the MSSM effective potential. 
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5 Conclusions 



In this paper we presented explicit and general results for the 0(ata s + o;|) two-loop corrections 
to the minimization conditions of the MSSM effective potential, which translate into corrections 
to the DR-renormalized parameters /i and m A . We discussed the numerical impact of our 
corrections in some representative scenarios of gravity-mediated SUSY breaking, and we found 
that the inclusion of the 0(ata s + ctf) corrections significantly improves the renormalization 
scale dependence of the results. Due to partial cancellations between the 0(ata s ) and 0(a^) 
corrections, including only the former may lead to inaccurate results. Our corrections are also 
required for consistency in the 0(a t a s + a|) two-loop computations of the MSSM Higgs masses, 
if the parameters entering the tree-level Higgs mass matrix are computed via renormalization 
group evolution from a set of high energy boundary conditions. 

A complete study of the electroweak symmetry breaking at the two-loop level would require 
also the knowledge of the corrections that are neglected in our gaugeless limit, among which 
the most relevant are controlled by the bottom Yukawa coupling and the electroweak gauge 
couplings. Concerning the corrections controlled by ctf,, they can be numerically non-negligible 
only for large values of tan/3. As discussed at the end of section^ the formulae for the 0(ai,a s ) 
corrections can be obtained by performing simple substitutions in their 0{ata 3 ) counterparts. 
On the other hand, the 0(a% + ottot^) corrections, which in some cases might be as relevant as 
the 0{abOi s ) ones, cannot be obtained in a straightforward way from the presently computed 
corrections and would require further work. 

In Ref. JZj a complete two-loop computation of the MSSM effective potential is presented, 
including also the terms controlled by the electroweak gauge couplings that are neglected in our 
analysis. The inclusion of such terms improves further the scale dependence of the parameters 
/i and m A , that in Ref. |17j are determined through a numerical minimization of the effective 
potential. However, the explicit analytical formulae for the two-loop tadpoles Si and S2, which 
are usually needed for practical applications, would be quite involved in the general case, and 
have not been presented so far. 

In conclusion, our work should lead to a more precise and reliable determination of the MSSM 
parameters at the weak scale, once the boundary conditions are provided at some larger scale 
according to the underlying theory of SUSY breaking. The fact that, among its many attractive 
features, the MSSM provides a natural mechanism for breaking radiatively the electroweak sym- 
metry, with the heavy top quark mass nicely falling in the required range 7 , seems to indicate the 
MSSM as the most viable theory for physics at the weak scale. However, only the forthcoming 
experimental results from the Tevatron and the LHC will tell us if this is indeed the case. 



7 As an example, we find that in the "typical" mSUGRA (SPSla) scenario the acceptable range for the top 
quark mass is 80 GeV < M t < 215 GeV. 
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Appendix A: Two— loop integrals 



We give here explicit expressions for the momentum integrals that appear in the two-loop part 
of the effective potential. The basic integrals in d = 4 — 2e dimensions are: 

,2(4-d) 



1 



(16vr 2 ) 2 
1 



J(x,y) 



(16 7T 



2^2 



I(x,y,z) 



(27T) 

2(4-d) 



2d 



d d p d d q 



(27T) 



2d 



[p 2 - x] [q 2 - y] ' 
d d p d d q 



[p 2 - x][q 2 - y][(p - q) 2 - z] 



(Al) 
(A2) 



Following Ref. PJJ], the functions J(x,y) and I(x,y,z) defined in Eqs. (|Alj) - (|A2j) are: 
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\ A(x, y, z) 

— — (x lnx + y Iny + z Inz) (lnx + Iny + Inz — 6) ^ — - — <£(x, y, z) . (A4) 

In the above formulae, lnx stands for ln(x/Q 2 ), where Q 2 = 47r^ 2 e _7B is the renormalization 
scale (je is the Euler constant). The functions A and <I> are respectively: 



A(x,y, z) = x 2 + y 2 + z 2 - 2 (xy + xz + yz) 



(A5) 



$(x,y,z) 



2 lnx_|_ lnx_ — Inu Inv — 2 I Li2(x+) + Li2(x_) I + 



7T 



(A6) 



where Li2(z) = — Jq dt [ln(l — t)/t] is the dilogarithm function and the auxiliary (complex) vari- 
ables are: 



x y 
u = — , v = — 



A = \J (1 — u — v) 2 — Auv , x± = - [1 ± (u — v) — A] 



(A7) 

The definition (|A6|) is valid for the case x/z < 1 and y /z < 1. The other branches of <3? can be 
obtained using the symmetry properties: 



$ (x, V, z ) = $ (y ; z> z) , x $ (x, y, z) = z $ (z, y, x) . 



(A8) 
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Finally, the following recursive relation for the derivatives of <3? proved very useful 8 for obtaining 
compact analytical results: 



. . . d $>(x, y,z) z 
y, z) = (y + z-x) ${x, y,z) + - 

OX X 



(y — z) In — h x ( In — + In — 



(A9) 



The derivatives of $ with respect to y and z can be obtained from the above equation with the 
help of the symmetry properties of Eq. (jA8|) . 



Appendix B: Explicit formulae for the 0(a t a s ) corrections 

We present here explicit expressions for the two-loop part of the functions F and G, giving rise 
to (D(ata s ) corrections to Si and £2. These formulae are valid when the parameters entering the 
one-loop parts of F and G are expressed in the DR scheme and computed at the renormalization 
scale Q. In units of g 2 Cf N c /(16tt 2 ) 2 , where Cf = 4/3 and N c = 3 are color factors, F 2e and 
G 21 read: 
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3 We thank G. Degrassi for explanations on how to derive Eq. ljA9|l . 
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The functions A(x,y, z) and <J>(x,y, z) appearing in F 2£ and G 2£ are defined in Eqs. (|A5|I - (|A6|) . 
The parameter S2Q is defined in Eq. Q13|). 
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tl,2 




t k,2 h,2 t 

(a) (b) (c) (d) 

Figure 1: Feynman diagrams that contribute to the two-loop effective potential and affect the 
0(ata s ) corrections to the electroweak symmetry breaking conditions. The diagrams relevant to 
the 0(aba s ) corrections can be obtained with the replacement t — > b. 




Figure 2: The classes of Feynman diagrams that contribute to the two-loop effective potential 
and affect the O(a^) corrections to the electroweak symmetry breaking conditions [q = (t,b), 
V = (H,h,G,A,H ± ,G ± ),h = (h°,h ± ), q = (tuhM)}. 
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"Typical" mSUGRA point 
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Figure 3: The DR parameters \fi(M z )\ (upper plot) and m A (M z ) (lower plot) as a function of 
the scale Q m i n at which the minimization conditions of the effective potential are imposed. The 
input parameters of the mSUGRA scenario are chosen as in the Snowmass Point SPS la |19| . 
The meaning of the different curves is shown in the caption and explained in the text. The mass 
of the top quark is taken to be 174.3 GeV. 
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"Coannihilation" mSUGRA point 
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Figure 4: Same as Fig. El for the Snowmass Point SPS 3. 
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"Large tanf3" mSUGRA point 
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Figure 5: Same as Fig. [21 for the Snowmass Point SPS 4. 
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'Light Stop" mSUGRA point 
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Figure 6: Same as Fig. El for the Snowmass Point SPS 5. 
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Figure 7: Same as Fig. El for the "Large mo" mSUGRA point proposed in section 21 
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